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Abstract 



Underactuated mechanisms provide low cost automation and can overcome actuator fail- 
ures. These mechanisms are more suitable for space applications mainly because of their less 
weight and lower power consumption. Typical examples of useful underactuated mechanisms 
in space would be large space structures and robot manipulators. Such mechanisms are how- 
ever difficult to control because of the fewer number of actuators in the system. In this paper 
we formulate the dynamics of an underactuated mechanism using Hamilton’s canonical equa- 
tions. Next, we develop a theorem that provides us with some necessary and some sufficient 
conditions for the asymptotic stability of autonomous systems. This theorem is more powerful 
than LaSalle’s theorem when higher order derivatives of the Liapunov function can be easily 
computed. Finally, we use a Liapunov function approach to develop a control strategy that will 
stabilize an underactuated mechanism in space to an equilibrium manifold. The effectiveness 
of such control is verified using our asymptotic stability theorem. 

1. Introduction 

Structures in space are mostly required for high precision tasks, like in the case of the 
orbiting interferometer telescope, or in the case of the space station that needs to point its 
antenna in a specific direction. These space structures are made up of trusses that are designed 
to have a light weight. The motivation behind this is to minimize the payload of the rocket 
that sends it in orbit. Light weight members have lower structural rigidity. Therefore while 
designing structural elements for space, much attention is paid to the geometric shape of the 
members so as to maximize their structural rigidity. Nevertheless, these trusses still posses 
a significant amount of structural flexibility. This is a serious disadvantage for large space 
structures because they easily pick up vibrations due to their flexibility. Vibrations could be 
thermally induced by differential heating of the structure or could be induced by differential 
gravitational forces. In the case of structures like the space station where robots are expected to 
perform routine tasks, vibrations could be easily induced through dynamic interaction between 
the robot and the structure. A free-flying multidegree of freedom system in space is a non- 
holonomic system (Nakamura and Mukherjee, 1990, 1991). Such systems have a noninvolutive 
property and they will experience a change in orientation under periodic motion (Kane, Head- 
rick and Yatteau, 1972; Vafa, 1987; Vafa and Dubowsky, 1987). Naturally, space structures 
will disorient themselves with time if vibrations persist. Though piezoelectric actuators may 
be used to damp out the vibrations in space structures, the system performs oscillations and 
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undergoes an undesirable change in orientation over a prolonged period of time. 

We now consider the prospect of replacing a single large flexible space structure with a chain 
of concatenated Light weight members. Each of the members of the chain can be considered 
rigid due to their smaller dimension and can be assumed to be concatenated with revolute 
joints. In order to achieve control over the system, we intend to use motors at some of these 
joints (instead of piezoelectric actuators). The other joints would be left unactuated. We intend 
to control such an underactuated mechanism such that it would be possible to configure the 
system in any desired way. If such control can be established, underactuated mechanisms would 
provide a meaningful alternative to large flexible space structures. 

Though the control of underactuated systems pose difficulties, in general they have a num- 
ber of advantages. Underactuated manipulators have lower power consumption and also weigh 
less. Therefore these manipulators will be very suitable for space applications. Besides space, 
underactuated systems w ill find applications for low cost automation, hyper-redundant manip- 
ulators, and manipulators with actuator failures. The range of tasks that can be performed by 
underactuated manipulators are however limited since these systems are usually incapable of 
exerting forces. This limitation can be overcome by the use of brakes at the unactuated joints. 
These brakes need not be used to stop the motion of the unactuated joints. Instead they may be 
used as clamps to maintain a fixed configuration of the unactuated joints over certain periods 
of time. These brakes would then enable the manipulator to perform tasks like force control. 
They would also allow the manipulators to behave as reconfigurable actuated systems. In the 
absence of brakes, underactuated manipulators may be used with proper control to pick and 
place objects, and to perform non-contact tasks like spray painting, arc welding, etc. 

Underactuated terrestrial robot manipulators were studied by Arai and Tachi (1990, 1991). 
In their studies, they assumed that the unactuated joints had brakes that could be used to 
stop the motion of the unactuated joints instantaneously. This simplification was used to 
eliminate the coupling between the actuated and unactuated links as and when desired. In 
1991, Arai and Tachi proposed a PID control law to control the trajectory of the actuated 
joints only. They verified the effectiveness of their control law through experiments on a 2DOF 
manipulator with one passive joint. Jain and Rodriguez (1991) studied the kinematics and 
dynamics of underactuated manipulators. They adopted the spatial operator algebra to develop 
an algorithm for the inverse dynamics. Papadopoulos and Dubowsky (1991) proposed the failure 
recovery control of space robotic systems. They showed in their formulation that it may be 
possible to control the joint whose actuator has failed when there exists a dynamical coupling 
between this joint and a joint whose actuator is functioning properly. Furthermore, in order 
that the passive joint can be controlled, the system inertia has to be invariant with respect to 
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the passive joint. 

In this paper we model our underactuated space mechanism with an open link tree structure 
consisting of (m + n) rigid links mounted on a free-flying space vehicle, n links of this mechanism 
are actuated while the rest m are unactuated. In the next section we use Hamilton's canonical 
equations (Goldstein, 1980; Nijmeijer and van der Schaft, 1990) to formulate the dynamics of 
this system. In section 4 we state and prove a theorem that provides us with some necessary and 
some sufficient conditions for the asymptotic stability of autonomous systems. This theorem is 
more powerful than LaSalle's theorem (LaSalle and Lefschetz, 1961) and gives us a systematic 
way to sort out the maximum invariant set from the set where the derivative of the Liapunov 
function (Liapunov, 1892) vanishes. Finally we use this theorem in section 5 to develop a 
control strategy for the stabilization of our underactuated mechanism in space to an equilibrium 
manifold 



2. Dynamics of free-flying underactuated systems 
- A Hamiltonian formulation 

In this section we formulate the dynamical equations of free-floating underactuated multi- 
body systems in space. We assume without any loss of generality that the system is of the form 
of a manipulator mounted on a space vehicle, as shown in Fig.l. We assume that the manipu- 
lator has a total of (m 4- n) joints, only n of which are actuated. The generalized coordinates 
of the system consist of q\ E R € representing the position and orientation of the space vehicle, 
<?2 £ representing the unactuated joint variables, and q 3 E R n representing the actuated 
joint variables. Due to the absence of gravitational potential energy in space, the Lagrangian 
Lo is equivalent to the kinetic energy of the system, and is given as 

Lo(q q) = \q T Mq, q = (qj qj qj ) T £ /? (6+m+n > ( 1 ) 

where, A 1 E /?( 6 + m + n ) x ( 6 + m + n ) is the inertia matrix of the system. It is a function of the joint 
variables q 2 and g 3 , but not a function of the vehicle variables q j. This is true because the 
kinetic energy' of the system is independent of the position and orientation of the space vehicle. 
Consequently, the Lagrangian is not a function of q\ and therefore the dynamics of the system 
is represented by the following vector equations: 



1 ( dL A 

dt \dqj 



= 0 



d_ 

dt 




( 2 ) 

(3) 
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where, T 6 R n represents the vector of the joint torques at the actuated joints. The right hand 
side of Eq.(2) is zero because we do not use the reaction jets or momentum wheels of the space 
vehicle. By refraining from using reaction jets, we can minimize the usage of jet fuel on board 
the spacecraft which is limited in quantity, and therefore we can maximize the useful lifespan 
of the system. 

When the Lagrangian is not a function of a set of generalized coordinates, like q : in 
our case, w'e call these coordinates cyclic or ignorable coordinates. In the presence of cyclic 
coordinates, some physical quantity of the system is conserved. In our case, the linear and 
angular momentum of the whole system is conserved. This conservation law is expressed by 
Eq.(2) and can be simplified to the form 



Miq = c (5) 

where A 1\ 6 /? 6x ( 6 + m + n includes the top six rows of the inertia matrix in Eq.(l), and C € 
R 6 represents the initial linear and angular momentum of the system. The above equation 
represents six velocity constraints on the motion of the system; three of these are holonomic 
while the other three are nonholonomic (Nakamura and Mukherjee, 1990, 1991). On the other 
hand, Eq.(3) represents m nonintegrable constraints that include second order derivatives of the 
generalized coordinates, and are therefore second order nonholonomic constraints. The degrees 
of freedom of our system are n, and is equal to the dimension of the control variable T. 

We now* use the transformation 

L(Q-Q,t) = Lo(q.q) + qjr ( 6 ) 

to define the input dependent Lagrangian function L(q,q,T ) (Nijmeijer and van der Schaft, 
1990). Under this transformation, we have the following relations 




By substituting Eqs.(6) and (7) into Eqs.(2), (3), and (4), we obtain the following homogeneous 
dynamical equations 



(9) 



d_ (dL_\ _ / 01 \ 

dt\dqj UgJ 

d_ ( 8L\ _ / 01 \ 
dt \dq 3 ) \dq 3 ) ~ 



( 10 ) 



We define the generalized momentum p E /?( 6 + m + n ) corresponding to the generalized 
coordinates q by the relation 






(ii) 



We next define the input dependent Hamiltonian function //(g,p, t) for the system, with 
the help of a Legendre transformation, as follows 



H(q.p.r) = p 1 q - L{q q.T 



( 12 ) 



Using the Legendre transformation in Eq.(12), the homogeneous dynamical equations given 
by Eqs.(8), (9). and (10) can be simply represented by Hamilton’s canonical equations 



(dH\ T fdH\ T 

Additionally, by substituting Eq.(6) in Eq.(12) we get the relation 



(13) 



H(q.p.r) = H 0 (q p)- qjr, H 0 (q,p) = p T q - L 0 (q.q) (14) 



which on differentiation yields 



H{q.p r) = Ho(q.p) - q 3 r - qj r 

° '• + (j?) T = H °- q = T -tfr (15) 

T 

By substituting the relation (dH/dr) = and the canonical expressions of Eq.(13) in the 
above equation, we finally get 



Ho (q,p) -qlr 



(16) 



We now go back to Eq. ( 14) for the definition of the Hamiltonian function Hq. Using 
Eqs.(l), (6), and (11) we can show that 



P = 




(17) 
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Therefore, from the definition of II o, 



H 0 = q T Mq - L 0 = ^q T Mq = L 0 (18) 

In other words, the Hamiltonian function Ho of our system represents the kinetic energy or 
equivalently the total internal energy of the system. Although Ho is equivalent to L 0 , it is a 
function of q and p only, and therefore the correct expression for Ho would be 

Ho = \p T M~ l p (19) 

which was obtained by substituting the relation q = A l~ l p from Eq.(17) into Eq.(18). 

3. Issues of stability and controllability 

In this section we first consider the stability of our nonlinear system from a linearization of 
the dynamics in the neighborhood of an equilibrium point. From Eqs.(13) and (14), our affine 
nonlinear system can be expressed by the form 

x = f(x) -f Br (20) 

x =(l)’ f = (dH 0 /dp dHo/dqf , B ^ J (21) 

where X 6 /? 2 ( 6 + m +") Is the state vector. / € R 2 ( 6 + m + n ),B e R 2{ and E n repre- 
sents the identity matrix of size n. Therefore a linearization of Eq.(20) around the equilibrium 
point (qo,Po-T 0 ) = (0.0,0) gives 



x = Ax -)- Br (22) 

(d 2 H /dpdq) M~ l \ r,2(6+m+n)x2(6+m + n) 

- (d 2 H/dq 2 ) - ( d 2 H/dqdp ) J 

From the definition of the matrix A it is clear that 

2(6+m+n) 

^ (dii/dxi) = tr(A) = 0 (23) 

* 

Equation (23) is the mathematical statement of Liouville’s Theorem (Goldstein, 1980). The 
above equation implies that the linearized system has as many eigenvalues in the open left half 
plane as those in the open right half plane. Therefore in the absence of the control vector T, 
we can conclude that the actual system is not exponentially stable. 
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The simplest approach to study the controllability of a nonlinear system as in Eq (20) is 
to consider its linearization. If the linearized system is found to be controllable, the nonlinear 
system is controllable in the neighborhood of the equilibrium point. However the lineariza- 
tion approach is often unsatisfactory. In the process of linearization the nonlinear system 
may loose much of its structure. Therefore a nonlinear system may be controllable though 
its linearization may not. In our case, it can be easily verified that the rank of the matrix 
( A — s.Z^2(m-f n+6) B ), where E t is the identity matrix of size i, is at most 2 (m + n) + 6 . 
Therefore the linearization of our system is not completely controllable. 

The controllability of a number of simple nonholonomic systems like the rolling contact 
(Li and Canny, 1990) and the single and multibody car systems (Laumond, 1987) have been 
individually studied by constructing the control Lie algebra. The control Lie algebra is defined 
as the smallest involutive distribution containing the span of the vector fields of the system and 
closed under Lie bracket operations. For these systems the local controllability was ascertained 
by showing that the rank of the control Lie algebra is equal to the dimension of the state space. 
It should be emphasized that unlike most of these nonholonomic systems, our system has a drift 
term (/ in Eq.(20)) due to the formulation of the problem at a dynamical level. Therefore the 
analysis based on the control Lie algebra cannot be performed on our system. 

In general our system may be asymptotically stabilizable by means of a linear or a non- 
linear feedback. However, Brockett (1983) has established some necessary conditions for the 
existence of smooth (infinitely continuously differentiable) stabilizing feedback laws for the gen- 
eral nonlinear system 



x = f(x,u), x € R* , u € R M , f(x t , 0) = 0 (24) 

with /(•,•) continuously differentiable in the neighborhood of the equilibrium point (:c f . 0 ). 
One of the three conditions require the mapping 

7 : R s x R si — R s defined by 7 : (x y u) f(x,u) 

to be onto an open set containing the origin, w’here X — 0. When q = 0, (dHo/dq) = 0 since 
Ho is quadratic in < 7 , and therefore from Eqs. (20) and (21) we arrive at 

p= (A) r (25) 

This clearly implies that the mapping 7 is not onto an open set containing the origin. Hence 
there does not exist a smooth feedback law' that can stabilize the system to an equilibrium point. 
This fact should how’ever not perturb us for we can always consider the problem of stabilizing 



the system to an equilibrium manifold, or stabilizing the system to an equilibrium point via 
a non-smooth feedback. In this paper we consider only the problem of stabilizing the system 
to an equilibrium manifold. In our next paper we shall address the problem of stabilizing our 
system to an equilibrium point via a non-smooth feedback. 

4. Theorem on Asymptotic Stability 

The Liapunov stability theorems provide sufficient conditions for proving the asymptotic 
stability of dynamical systems. For autonomous systems these theorems are easy to apply 
when we can show’ that the derivative of the Liapunov function is negative definite. When the 
derivative of the Liapunov function is negative semidefinite, we often face problems. In such 
situations it may be possible to conclude the asymptotic stability of the system using LaSalle's 
theorem provided we can show that the maximum invariant set contains only the equilibrium 
point. It is always possible and easy to identify the set of points w’here the derivative of the 
Liapunov function vanishes but the maximum invariant set is only a subset of this set. The 
main challenge of LaSalle’s theorem is therefore to sort out the maximum invariant set. More 
importantly, LaSalle’s theorem is inapplicable to nonautonomous systems. In the event where 
the derivative of the Liapunov function vanishes, there exists no readily applicable result for 
proving the asymptotic stability of nonautonomous systems. 

In this section we develop sufficient conditions for proving the asymptotic stability of 
autonomous systems when the first derivative of the Liapunov function is negative semidefinite. 
These sufficient conditions involve higher order derivatives of the Liapunov function that contain 
information of the higher order dynamics of the system. Consequently, it becomes easier to 
identify the maximum invariant set. In this section we also provide some necessary conditions 
for the asymptotic stability of autonomous and nonautonomous systems. Before stating our 
asymptotic stability theorem, we state the following Lemmas. 

Lemma 1 . A real function /(f) 6 C 2 defined in (a, 6) is concave iff /"(f) < 0, Vf 6 (a, 6). 
Lemma 2. Let /(f) be a nonpositive function such that /(fo) = 0 and /(f) < 0 for some values 
of f. If the function /(f) is analytic, then /(f) is concave is some open neighborhood of f 0 . 

The proofs of the two Lemmas stated above have been provided in the Appendix for reference. 
Using these two lemmas we can conclude that if /(f) is a smooth nonpositive function and 
/(f o) = 0, then /'(fo) = 0 because /(f) is locally maximum at f 0 , and /"(f) < 0 in some open 
neighborhood of f 0 . If /"(fo) = 0 also, then we can apply our lemmas to /"(f). In such a case 
/'"(fo) = 0, and /""(f) < 0 in some open neighborhood of f 0 . Our lemmas can therefore be 
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applied recursively. When some even derivative of f(t) vanishes at / 0 , the next higher derivative 
which is an odd derivative also vanishes at t o, and the second next derivative is nonpositive in 
some open neighborhood of to- 

Let us now consider the nonautonomous system 

x = f(t,x) (26) 

where f : R+ x D — ■ /? n is a smooth vector field on R+ x D, D C R n is a neighborhood of 
the origin X = 0. Let X — 0 be an equilibrium point for the system described by Eq.(26). We 
then have 



/(/, 0) = 0. V/ > 0 (27) 

We next state and prove our theorem on asymptotic stability. 

Theorem 1. Let : R+ x D — R+ be locally positive definite and smooth on R+ x D , 

such that 

i’(tz) = BV/dt + (dV/dx) f(t,x) (28) 

is locally negative semidefinite. Then whenever an odd derivative of \ ' vanishes, the next 
derivative necessarily vanishes and the second next derivative is necessarily negative semidefi- 
nite. Furthermore, a sufficient condition for an autonomous system to be asymptotically stable 
is that there exists a positive integer k such that 

f r ( 2 i +'>(z) < o vx : i'(i) = o f2q 

\r<‘>(x) = 0 for i = 2,3,---, 2k K 1 

where r (,, (x) denotes the (*) — th time derivative of l' with respect to time. 

Proof: The necessary conditions of this theorem can be proven very easily with the help of 
Lemmas 1 and 2. 

To prove that Eq.(29) provides sufficient conditions for asymptotic stability, we first realize 
that X = 0 is stable by standard argument since V is locally positive definite and V < 0. 

Next, since V is bounded from below by zero and V is nonincreasing (1* < 0), V — * a, 
a > 0, as t — ► oo. 

Since V is smooth, V is uniformly continuous. Hence when V — ► a, V — ► 0 as / — ► oo, by 
Barbalat’s lemma (Slotine and Li, 1991). 

Since V is locally positive definite, — 0 => X — ► 0 as f — oo. Therefore if we can 

show that q = 0 w r e can conclude asymptotic stability. We prove a = 0 by contradiction. Since 
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V — > a ± 0 and V is locally positive definite, 3 an open neighborhood A' of x = 0 such that 
the trajectory of x(t) lies outside A r V/ > T, and for some T > 0. 

Let 5 = ja: : V’(x) = oj. Since ar(f) converges to 5 but lies outside A T for large the set 
W = S - A r is nonempty and is the limit set for x(f). Then let U be an open neighborhood of 
\V whose closure U c does not contain X — 0. Now let us denote 

—7 = max V( 2k + l \x) (30) 

x€U c 

Then —7 < 0. Since x(t) — < ► W as t — ► 00 , 3 T\ such that x(t) £ U c > T\. Now integrating 
V /r ( 2t + i)(t) with respect to time to get V, we have 



V(t) - V(T \ ) = / ... / 
7r, v/t, 



(2* + l) 



</ /'- 

Ti v/Ti 




(« + .,! (31 » 

Hence V(0 < V(7i) - 7 (< - 7i) 2t+1 /(2* + 1)!. Since V(Ti) < V(t = 0), V(t) oo as 

t — > oo. This contradicts the fact that V > 0. Hence a = 0 and that implies that the system is 
asymptotically stable. 



5. Stabilization to an equilibrium manifold 

The state variables of a free-flying underactuated mechanism in space was shown in Eq.(21) 
as x = ( q T p T ) T , where q = ( qj q \ qj ) T £ /? 6+m+n denote the generalized coordi- 
nates and p £ J R 6+m+n denote the generalized momentum corresponding to these generalized 
coordinates. We intend to control the system in such a way that p — ► 0 and q$ — * q$d at the 
final point of time. q$d denotes the desired configuration of the actuated joints of the system. 
If such a control can be established, the underactuated system would come to a complete rest 
and the actuated joints of the system w ? ould converge to their desired values simultaneously. 
We therefore define a Liapunov function (Liapunov, 1892) v as 

V = Ho + 1 Aql A q 3 , A q 3 = (q 3 j - q 3 ) (32) 

where Hq is the Hamiltonian of the system defined by Eqs.(18) and (19). Since the Hamiltonian 
Hq represents the total kinetic energy of the system. Hq = 0 is attained only when p = 0 or 
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alternatively q = 0. If we now define an equilibrium manifold M c — {x q 3 = q 3d P = 0}. 
then the Liapunov function r defined by Eq.(32) is zero only on the equilibrium manifold and 
positive everywhere else. The derivative of t is next computed as 



t = H 0 - A ql q 3 

= qlT-AqJq 3 = ql(T-Aq 3 ) (33) 

w’here Ho = qj T was substituted from Eq.(16). We now choose T in Eq.(33) as 

r = Ag 3 - 0Q3 (34) 

where 0 is a positive scalar quantity. Substitution of Eq.(34) in Eq.(33) yields 

i' = -0\\ 93 || 2 (35) 

Clearly, r is negative serrudefinite and is equal to zero w ? hen q 3 = 0. At this point LaSalle's 
theorem (LaSalle and Lefschetz, 1961) could be used to conclude the asymptotic stability of the 
system to the equilibrium manifold provided we could show’ that q 3 = 0 is attained only when 
<73 = Q 3 d and p = 0. Since LaSalle's theorem does not provide us with any systematic way to 
sort out the maximum invariant set from the set of all X : v = 0, w’e refer to our theorem that 
was stated and proved in the earlier section. 

By computing the second and the third derivatives of the Liapunov function r from Eq.(35) 
w'e can show’ that w’hen 1 = 0 or equivalently q 3 — 0, r = v ( 2) = 0 and = — 2 0\\ q 3 || 2 < 0. 
Additionally if C 31 = 0 then (73 = 0. Then we can show by computing the higher order 
derivatives of v that t = 0 and r (5) = -6/?|| q 3 3) || < 0, w’here g3 3> is the third derivative of 
<73 wdth respect to time. In other words whenever an odd derivative of the Liapunov function r 
vanishes, the next derivative also vanishes and the second next derivative is found to be negative 
semidefinite. This is in complete agreement wdth Lemmas 1 and 2. Furthermore this satisfies 
the necessary conditions of our asymptotic stability theorem. 

From the above discussion it follow’s that the choice of the control vector T in Eq.(34) 
results in 



f( 2l+| ) = -AH q< 3 k+i) || 2 for > 0, andfor k = 1, 2, • ■ • , oo (36) 

when t>(‘) = 0 for i — 1,2,---, 2 k. Therefore when v = 0 or equivalently q 3 = 0, if q 3 i + l ’ / 0 
for some positive integer k , then the sufficient conditions of our theorem given by Eq.(29) are 
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satisfied and we can conclude asymptotic stability of our system to the equilibrium manifold 
Me- 



6. Conclusion 

We have discussed in our paper the dynamics and control of underactuated mechanisms 
in space. The dynamics of the system was formulated using Hamilton’s canonical equations. 
To prove the stability of our system we have developed a general asymptotic stability theorem. 
It is an elaboration of LaSalle’s theorem and it provides us with a systematic way to sort out 
the maximum invariant set from the set where the dervative of the Liapunov function vanishes. 
Similar to LaSalle’s theorem, the limitation of our theorem is that it is applicable only to 
autonomous systems. Using a Liapunov function approach we have developed in this paper a 
control strategy that brings an underactuated mechanism to rest and converges the actuated 
joints to their desired configuration simultaneously. We show that our control law is effective 
provided all the derivatives of the actuated joint velocities are not zero simultaneously. Hence 
the observability of the system, when we take the velocity of the actuated joints as the output, 
is going to play an important role in the stabilization. This part is going to be w T orked out 
using a geometric nonlinear control approach and will be appended to the paper by the time it 
is due. We will also include examples of situations where observability is lost, and provide the 
physical meaning of such situations. 
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7. Appendix 



Lemma 1. A real function /(/) £ C 2 defined in (a, 6) is concave iff /(f)" < 0, Vx £ (a. 6). 



Proof: 

(a) Necessity 

Let x £ (a, 6). Then for h small enough, x — h, x + h £ (a, 6). From the definition of concavity 
(Rudin), f(x) > i (f{x - h) + f(x + h)). Therefore, since / £ C 2 , 



f'(x) = limh^o 



/(x-/t) + /(x + fi)-2/(x) 

h 2 



< 0 



(.4-1) 



(b) Sufficiency 

Let x, y £ (a.b). and x < y. For A £ [0, 1], and t = Ax + (1 — A )y, the first order Taylor’s series 
approximation of /(x) and f(y) are respectively 



/(*) = f(t) + f'(t)(x-t) + f"(ti)(x-t) 2 , *,€[*,*] (A- 2) 

f(y) = m + f'(t)(y-t) + n^)(y-t) 2 , 6€[f,y] (A- 3) 

Therefore it follows that 



A/(x) + (1 - A)/(y) = f(t) + A /"(6)(x - f) 2 + (1 - A)/"(fe)(y - f) 2 

</(f) since /"K,) < 0, /"(6)<0 (A- 4) 



Therefore the function is concave by definition. 



Lemma 2. Let /(f) be a nonpositive function such that f(to) — 0 and f(t) < 0 for some values 
of t. If the function f(t) is analytic, then f(t) is concave is some open neighborhood of f 0 - 
Proof: Since the function f(t) is analytic, all derivatives of the function exist and the function 
can be expanded using Taylor’s series as 



/«) = E - «.)” M - 12) 

n = 0 

Let us next assume that our function f(t) is not concave in any open neighborhood of t 0 . This 
implies from Lemma 1 that the condition f"(t) < 0 does not hold good in any open neigh- 
borhood of t 0 . Therefore either /"(f) > 0, or /"(f) changes sign in every open neighborhood 



14 



or 1.. If no > 0 every open neighborhood of then we can chow from the cotollatv of 
Lemma 1 that /,<) » convex everywhere. Th, s is no, true because /(,) is noupostttve and'has 

a maximum value at , - The other possibility is that /"(,) changes sign in every open 

neighborhood of t 0 . Then f {n) (t) for n - 2 1 . r*- „ 

j W iom_^j,... iOC changes sign in every open neighborhood 

of < 0 - This implies that f (n) (t n ) = 0 for n - 2 1 ~ ajjv „ 

n ~ Additionally, since f(t) is nonpositive 

and /(«„) = 0, /(,) achieves a local maximum a. Therefore /'«„) = 0 Substituting these 
results Eq.(A-12), have /(,) = 0. Thus cannot be „ue because /(,) is strictly negative fo, 

some values of I. We have therefore proved by contradiction that /(,) is concave some open 
neighborhood of t 0 . 
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Fig.l. A free-flying under-actuated mechanism in space 
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